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Abstract 

Unitarity constraints on the scalar parameters both for the Standard Model 
and the general Two Higgs Doublet Model (THDM) are examined. In the case 
of the THDM with an exact discrete symmetry transformation, we show that 
the mass of the lightest CP-even Higgs boson (Mh) and tan (5 are strongly 
correlated and consequently a strong lower bound can be put on Mh for large 
tan (3. It is also shown that the inclusion of the discrete symmetry breaking 
term relaxes the aforementioned correlation. 
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1 Introduction 



The Standard Model (SM) of electroweak interactions U is in complete agreement 
with all precision experimental data (LEP,Tevatron, SLD). What remains to be 
discovered is the standard Higgs boson (0°), which is the major goal of present 
and future searches at colliders. From the 95% CL upper limit obtained from the 
global Standard Model fit 0, it seems that the existence of a relatively light Higgs 
boson is favored (M^o < 210 GeV at 95% confidence level), which is a feature of 
Supersymmetric theories. If we believe in SM, the global Standard Model fit tell 
us that Higgs boson particle is going to be seen very soon. Direct searches at LEP 
show no clear evidence for the Higgs boson, and imply a lower limit tjih > 112.3 
GeV at the 95% CL || (Note that this is still a preliminary limit, the previous 
combined limit is m# > 107.9 GeV). Last September exciting new results from 
LEP experiments has postponed by one month the LEP program. At the end of 
this extra period, the data collected by the four experiments of LEP are compatible 
with a production of Higgs boson in association with Z boson but also compatible 
with other known process. The four experiment of LEP agreed that a further run 
with about 200 pb -1 per experiment at a center-of-mass energy of 208.2 GeV would 
be enable the four experiments to establish a 5a discovery Q. If there is no extra- 
extension for LEP run, we have to wait until either the Runll at Tevatron || may 
aid in the interpretation of the signal or the LHC, which is able to discover such a 
light Higgs, start up 0. 

Constraints on M^o can be obtained by making additional theoretical assump- 
tions. These include i) the triviality bound [0, ii) Vacuum stability requirement and 
iii) demanding that unitarity of the S-matrix is not violated. 

The triviality bound follows from the fact that the quartic coupling of a pure 
A$ 4 scalar theory increases as a function of the momentum scale, causing the theory 
to become non-perturbative at the so-called Landau pole. The A$ 4 theory will then 
have X(Q) = at some scale Q which renders the theory trivial. 

By requiring at some large scale A (where new physics is supposed to enter) 
that the quartic coupling A (A) satisfies the relation -w^y > (which means that 
A (A) remains non- vanishing and so the theory is non-trivial) one can obtain an 
upper bound on the Higgs boson mass. For example, if this large scale A approaches 
the grand unification scale 10 15 GeV one can find a strong approximate upper bound 
on the Higgs boson mass || M^o < 160 GeV [] 

The condition of vacuum stability requires that A remains positive at a large 
scale in order to have spontaneous symmetry breaking; if A becomes negative then 
the potential is unbounded from below and has no minimum. Including the 2 loop 
renormalisation group improved effective potential, the vacuum stability require- 
ment gives (if the Standard model is valid up to scale 10 16 GeV) f|, |§: M<f>° > 
130.5 + 2.1 (M t — 174) GeV. Finally, unitarity arguments lead to an upper bound on 
the Higgs boson mass |l(J, and this will be covered in more detail in section 2. 

In recent years there has been growing interest in the study of extended Higgs 

1 The above bound changes if we couple the theory to fermions and gauge bosons. 
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sectors with more than one doublet JTT|. The simplest extension of the MSM is the 
Two Higgs Doublet Model (THDM), which is formed by adding an extra complex 
SU(2)i ® U(l)y scalar doublet to the MSM Lagrangian. Motivations for such a 
structure include CP-violation in the Higgs sector, supersymmetry, and a possible 
solution to the cosmological domain wall problem jT^]. In particular, the Mini- 
mal Supersymmetric Standard Model (MSSM) |TTJ takes the form of a constrained 
THDM. The two most popular versions of THDM are classified as type I and type II, 
differing in how the Higgs bosons are coupled to the fermions although both versions 
possess identical particle spectra after electroweak symmetry breaking. From the 
eight degrees of freedom initially present in the two Higgs doublets, three correspond 
to masses of the longitudinal gauge bosons, leaving five degrees of freedom which 
manifest themselves as five physical Higgs particles: Two charged Higgs H , two 
CP-even H°, h° and one CP-odd A . Until now no Higgs boson has been discov- 
ered, and from the null searches one can derive direct and indirect bounds on their 
masses. The latest such limits are: for the charged Higgs boson the LEP combined 
result is m#± > 77.5 GeV 0, and for the neutral Higgs OPAL collaboration [[TJ| 
has made the following scan of the parameter space as: 1 < < 100, 5 < rriA < 2 
TeV, — 7r/2 < a < and 0.4 < tan/3 < 58 finding that the region ( 1 < m h < 44 
GeV and 12 < rriA < 56 GeV ) is excluded at 95% independent of a and tan/3. 

We note here that tree-level unitarity constraints for the THDM scalar poten- 
tial constrained by the exact discrete symmetry <3>j — > — were studied in W\ [15 



When deriving constraints from unitarity, 0] considered only seven elastic scatter- 
ing processes S1S2 — > S1S2 (where Si is a Higgs scalar) while [T5|] considered a larger 
scattering (S) matrix. In [fT5| , upper bounds on the Higgs masses were derived, in 
particular Mh < 410 GeV for tan/3 = 1, with the bound becoming stronger as 
tan/3 increases. In section. 3, at the same footing as we will present for the MSM 
in section. 2, we improve the aforementioned THDM studies by including the full 
scalar S matrix which includes channels which were absent in |Tj|. In addition, we 
also show graphically the strong correlation between Mh and tan (3 both in the case 
where the discrete symmetry $3 — ► — $j is exact and also in the case where it is 
broken by dimension 2 terms: X^(^f^2)- 

The paper is organized as follows. In Section 2 we present the unitarity ap- 
proach we will be using and review the unitarity constraints in the framework of the 
MSM. Section 3 contains a short review of the THDM potential and the analytical 
unitarity constraints. In Section 4 we present our numerical results for the cases of 
A5 = and 7^ 0, while Section 5 contains our conclusions. The appendix contains 
pure scalar quartic interactions both for MSM and for THDM. 



2. Unitarity approach 

We will study the unitarity constraints by computing the scalar scattering processes, 
S1S2 — > S3S4, i.e. both elastic and inelastic channels. In terms of the partial wave 
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decomposition, the amplitude M. of a scattering S1S2 — > S3 S4 can be written as: 

00 

M(s,t,u) = 167r^(2Z + l)P(cos#)a,(s) (2.1) 
1=0 

Where s, t, u are the Mandelstam variables, ai(s) is the spin I partial wave and Pi 
are Legendre Polynomials. The differential cross section for S1S2 — > S3S4 is given 
by: 

— = —\M\ 2 
dtt 64tt 2 s 1 1 

Using the fact the the Legendre polynomials are orthogonal, the cross section be- 
comes: 

1 f)7T 

= E£ (2Z + l)h| 2 

s 

The optical theorem together with expression for the total cross section given above 
leads to the following unitarity constraint: 

^{aif + % ai ) 2 = \ai\ 2 = 9f(a,) for all I 

The above equation is nothing more than an equation for a circle in the plane 
(9ft(aj), 3(a/)) with radius ~ and center (0, |). It can be shown easily from the 
graphical representation of this circle that: 

|£(a,)|<2 for all I 

The partial wave aj(s) can be inversed from eq.^TT] and one finds: 

1 /-i 

a/(s) = / ci(cos6')P;(cos6')A<(s, t, u) 

If we limit ourselves to the J = s-wave amplitude ao(s), for vanishing external 
masses of Si^,3,a, flo(s) take the following form: 

aoOO = iJ-[Q + { WrAjS " J feW + cJiWJMl + 4)H2-2) 

107T S h 

Where Q is the four point vertex for SiS% — > S3S4 and are the trilinear vertex 
hSiSj, q = 1 or (resp c u = 1 or 0) for processes with or without t-channel (resp 
for process with or without u-channel). 

The first term in the bracket in the r.h.s of eq.2^ is the contribution of the quartic 
coupling Q = S1S2S3S4 to the amplitude Ai; the second term is the contribution 
of the s-channel diagram with exchange of particle h and the third term is the 
contribution of t and u channels. 



In very high energy collisions, it can be shown from eq j2.2| that the dominant 
contribution to the amplitude of the two-body scattering SiS 2 — > S3S4 is the one 
which is mediated by the quartic coupling. Those contributions mediated by tri- 
linear couplings are suppressed on dimensional grounds. Therefore the unitarity 
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constraint \a Q \ < 1/2 reduces to the following constraint on the quartic coupling, 
iQiSiS^S^Si)] < 87T. In what follows our attention will be focused on the quartic 
couplings. 

To constrain the scalar potential parameters one can demand that tree-level 
unitarity is preserved in a variety of scattering processes. This corresponds to the 
requirement that the J = partial waves (do) for scalar-scalar, gauge boson-gauge 
boson and gauge boson-scalar scattering satisfy |ao| < 1/2 in the high-energy limit. 
At very high-energy, the equivalence theorem |[TjJ states that the amplitude of a 
scattering process involving longitudinal gauge bosons V^'° may be approximated 
by the scalar amplitude in which gauge bosons are replaced by their corresponding 
Goldstone bosons G ±,Q . We conclude that unitarity constraints can be implemented 
by solely considering pure scalar scattering. 

Let us apply all these features to the non-coupled process W^W£ — > W^W^ 
in the MSM, where the subscript L denotes the longitudinal polarisation states. As 
stated above, at high energy the amplitude of W^W£ — > W^W£ is approximated 
by G + G~ — > G + G~ (where G is the charged goldstone associated with W ) whose 
dominate contribution comes the quartic coupling, G + G~G + G~ = M| /f 2 , where 
v is fixed by the electroweak scale as v 2 = (2\/2Gf)~ 1 - The unitarity constraint 
\Q\ < 87r gives the following upper bound on the Higgs boson's mass: 



Mj < 8nv' 



M 2 LQT 



4tt 



V2Gi 



(870) 2 GeV 2 



Where Mlqt is the bound deduced by Lee, Quigg and Thacker [TIJ . 
In fact the channel W^W^ considered above is coupled with the following chan- 
nels: ZlZl/\/2, 0°0°/\/2 and Z L (f)° (the factor a/2 accounts for identical particle 
statistics). Together with those four neutral initial states we consider also the 2 
charged channels W^(p and W^Z. Note that due to charge conservation the two 
charged channels are not coupled to the four charged channels. Taking into account 
all the above channels, the scattering amplitude is given by 6 x 6 matrix which 
is diagonal by block: 4x4 block for the 4 neutral channels and 2x2 block for 
2 charged channels. At high energies, the matrix elements are dominated by the 
quartic couplings and so the full matrix in the basis {W^W£ , ZlZl/\/2, </> </> /\/2, 
Z L (p ,W£(j> ,W£Z) takes the following form: 



M 2 



/ 


1 


V2 
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\) 



(2.3) 



This matrix can be read easily from the appendix A.l in which we have listed the 
pure scalar quartic interactions in the MSM. 

The requirement that the largest eigenvalues^ of ao respects the unitarity constraint 
2 In the case of this matrix the largest eigenvalues is 3/2 
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M 2 < -8ttv 2 = -M 2 LQT « (710) 2 GeV 2 



yields |I7||: 

2 2 2 

3 3 

We conclude that the inclusion of the complete set of scattering channels (charged 
and neutral ones) into the analysis, leads more or less to a stronger unitarity con- 
straint on the Higgs boson mass M^o < 710GeV. The unitarity requirement tells 
us that if the Higgs boson mass is above the unitarity constraint then the standard 
model becomes non-perturbative and some new physics should appear at some high 
scale to restore the unitarity of the theory. 



3. THDM scalar potential and Unitarity constraints 
3.1 THDM scalar potential 

The most general THDM scalar potential which is renormalizable, gauge invariant 
and CP invariant depends on ten parameters, but such a potential can still break CP 
spontaneously Jl8|]. In order to ensure that tree-level flavor changing neutral currents 
are eliminated, a discrete symmetry ($j — > — where is a scalar doublet) may 



be imposed on the lagrangian [I9[, which reduces the number of free parameters to 



6. The resulting potential was considered in |p0| , and is referred to as Va in |2l 
We shall be concerned with the potential described in which is equivalent to Va 
plus a term which breaks the discrete symmetry (parameterized by A5) and contains 
7 free parameters. Such a potential does not break CP spontaneously or explicitly 
20[],[22[ provided that all the parameters are real. 



It has been shown that the most general THDM scalar potential which is 
invariant under SU(2)l <S> U(1)y and conserves CP is given by: 

V(^, $ 2 ) = Axd^l 2 - v\) 2 + A 2 (|$ 2 | 2 - v 2 2 f + A 3 ((|$ i| 2 - vl) + (|$ 2 | 2 - vl) f + 

A 4 (|$i| 2 |$2| 2 - \®xH 2 ) + A 5 (&($f$ 2 ) - V1V2) 2 + A 6 [3($+$ 2 )] 2 (3.1) 

where $1 and $ 2 have weak hypercharge Y=l, V\ and v 2 are respectively the vacuum 
expectation values of $1 and $ 2 and the Aj's are real-valued parameters. 



1 hj+izj 



This potential violates the discrete symmetry $j — > — $j softly by the dimension 2 
term As^^^) and has the same general structure of the scalar potential of the 
MSSM. One can prove easily that for A5 = the exact symmetry $j — > — $j is 
recovered. 

After electroweak symmetry breaking, the W and Z gauge bosons acquire 
masses given by m 2 ^ = ^ 2 f 2 and m| = |(g 2 + g' 2 )v 2 , where g and g' are the 
SU (2) weak and U (l)y gauge couplings and v 2 = v\ + v\. The combination v\ + v\ is 
thus fixed by the electroweak scale through v \+v\ = (2v / 2Gi?) _1 , and we are left with 
7 free parameters in eq.( |3.1|) , namely the (Aj)i=i i .„ ) 6's and tan/3 = f 2 /fi. Meanwhile, 
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three of the eight degrees of freedom of the two Higgs doublets correspond to the 3 
Goldstone bosons (G , G°) and the remaining five become physical Higgs bosons: 
H°, h° (CP-even), A (CP-odd) and H . Their masses are obtained as usual by 
the shift $i — >• $j + Vi [|TT| and read 

1, 



m 2 A0 = X 6 v 2 ; m 2 H± = \ 4 v 2 and m 2 H0h0 = -[A + C ±y(A — C) 2 + 4B 2 ][3.2) 
where 

A = Avf(\x + A 3 ) + V2X5 , B = v x v % {AX z + A 5 ) and C = 4v% (A 2 + A 3 ) + v 2 \ 5 

The angle (3 given by tan/3 = 1)2/ Vi defines the mixing leading to the physical H ± 
and A states, while the mixing angle a associated to H°, h° physical states is given 



by (See JTT| for details.) 



2B A-C . , 

sin 2a = . , cos 2a = . (3.3) 

'(A-Cy + AB 2 J {A - C) 2 + AB 2 



It will be more suitable for the forthcoming discussion to trade the five parameters 
Ai, 2,4,5,6 for the 4 Higgs masses and the mixing angle a. From now on we will take the 
physical Higgs masses, m H o , m^o , m^o , m H ± , the mixing angles a, (3 and the coupling 
A5 as the 7 free parameters. 

It is then straightforward algebra to invert equations (ft.2|) through (|3~3|), and 
get the Aj's in terms of this new set of parameters. 

\ 9 2 2 \ 9 2 2 \ 9 1 s a c a / 2 2\ Ak. .x 

K = ^ mH± ' Ae = 2<;^ ' A3 = 8<^ (m "- m ^ ^ i 3 - 4) 

2 \ 

A i = o 2 2 i c l m H + s a m fe - r^^H-wifc)] - ^(-1 + tan 2 /?) (3.5) 
oc^TTi^y tan jj 4 



a 2^ 2 i s l m H + c l m l ~ s a c a tan (3(m 2 H - m 2 h )\ - ^(-1 + — 



3.2 Unitarity constraints 

To constrain the scalar potential parameters of the THDM one can demand that 
tree-level unitarity is preserved in a variety of scattering processes: scalar-scalar 
scattering, gauge boson-gauge boson scattering and scalar-gauge boson scattering. 
We will follow exactly the technique developed in section. 2 and therefore we limit 
ourselves to pure scalar scattering processes dominated by quartic interactions. 

In order to derive the unitarity constraints on the scalar masses we will adopt 
the technique introduced in It has been shown in previous works \2E\ that the 



quartic scalar vertices written in terms of physical fields H ± , G ± , h°, H°, A and G°, 
are very complicated functions of Aj, a and (3. However the quartic vertices (com- 
puted before electroweak symmetry breaking) written in terms of the non-physical 
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fields wf , hi and Zi (i= 1,2) are considerably simpler expressions. The crucial point 
of [15 is the fact that the S matrix expressed in terms of the physical fields (i.e. 



the mass eigenstate fields) can be transformed into an S matrix for the non-physical 
fields <f i , hi and Zi by making a unitarity transformation. The latter is relatively 
easy to compute from eq, 



3.1 



Therefore the full set of scalar scattering processes 
can be expressed as an S matrix 22 x 22 composed of 4 submatrices [ A^i(6 x 6), 
M 2 {Q x 6), ^3(6 x 6) and ^4(8 x 8)] which do not couple with each other due to 
charge conservation and CP-invariance. The entries are the quartic couplings which 
mediate the scattering processes. The pure scalar quartic interactions expressed in 
terms of the non-physical fields hf, hi and Zi are listed in the appendix A. 2. 

The first submatrix M\ corresponds to scattering whose initial and final states 
are one of the following: (wiW2,w 2 Wi , hiz 2 , h 2 Zi, Ziz 2 , hih 2 ). With the help of 
appendix A. 2, one can find that M\ takes the following form: 



/ 2A 3 + A 56 
A 



Mi 



/2 



A 



-*A 6 - 4 /2 



iX^/2 



56 

iX^/2 
-iX^/2 
A54/2 
A54/2 

— A; i A, 1 . 



V 

where Xf, 
following 5 distinct eigenvalues: 

e x = 2A 3 
e 2 

.A 



2A 3 + A+ /2 zA 6 " 4 /2 -iX^/2 

-iX^/2 2A 3 + A 6 A 5 - 6 /2 

iX^/2 XVJ2 2A 3 + A 6 

A54/2 

A54/2 



A54/2 A54/2 \ 

A54/2 A 54 /2 





2A 3 + A 5 Ase/2 

Ase/2 2A 3 + A 5 J 



r With the help of Mathematica we find that Mi has the 



A A5 + 5 A 
A 4 - — + -A 6 

2A 3 + A 4 — — + -A 6 



2A, 



A4 + 2 As 



/- — 2A3 + A4 + -A5 — -Xe 
fi — / 2 = 2A 3 + -A 5 + -A 6 



(3.7) 



The second submatrix M 2 corresponds to scattering with initial and final 
states one of the following: (wfwi , w 2 w 2 , ^p), where the \[2 

accounts for identical particle statistics. Again, with the help 01 appendix A. 2, one 
find that A^2 is given by: 



M, 



( 4A+ + 


2A 3 + A+ 12 


V2Xf 3 


X34/V2 


V2Xf 3 


~X 3i /V2 \ 


2A 3 + ^f 


4A 2 f 3 


X 3i /V2 


V2X+ 


~X 3i /V2 


V2X+ 


V2A+ 


WV2 


3A+ 3 


A35/2 




A 36 /2 


~X M /V2 


V2X+ 


A35/2 


3A 2 3 


A 36 /2 


\ + 

A 23 


V2X+ 


~X 3 a/V2 


Ai 


A 36 /2 


3A 13 


A35/2 


\ X34/V2 


V2X+ 


A 36 /2 


A23 


A35/2 


3A 2 3 / 



(3.8) 
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where A3i = 2 A3 + Aj. The matrix M2 possesses the following 6 distinct eigenvalues: 



by: 



a± = 3(Ai + A 2 + 2A 3 ) ± W9(A! - A 2 ) 2 + (4A 3 + A 4 + ^(A 5 + A 6 )) 2 (3.9) 



b ± = Ax + A 2 + 2A 3 ± J(Ai - A 2 ) 2 + i(-2A 4 + A 5 + A 6 ) 2 (3.10) 



c± = Ai + A 2 + 2A 3 ± J (Ax - A 2 ) 2 + ^(A 5 - A 6 ) 2 (3.11) 
The third submatrix Ai 3 corresponds to the basis: (hiZi, h^z-i) and is given 



" M 3=f^ 3 (3-12) 

\ 2 A 56 ZA 23 / 

The matrix .M3 possesses the eigenvalues d± and c±, with d-t = c±. All the above 



eigenvalues agree with those found in Ul5 |, up to a factor of 1/167T which we have 



factorised out. In our analysis we also include the two body scattering between the 
8 charged states: hiWi, h 2 Wi, ziwf, z 2 Wi, hi%u 2 , h 2 w 2 , Z\W 2 ', z 2 w 2 . The 8x8 
submatrix M.4 obtained from the above scattering processes is given by: 



M 4 = 



( 2A+ 














A54/2 





-zA 6 " 4 /2 





2A 3 + A 4 








1 \- 

2 A 54 





^2 A 64 











2A 1 f 3 








1 \~ 

2 A 64 





A54/2 











2A 3 + A 4 


— 2 A 64 





1 \- 

2 A 54 








1 \- 

2 A 54 





2 A 64 


2A 3 + A 4 











1 \- 

2 A 54 





2 A 64 








2A23 











2 A 64 





1 \- 

2 A 54 








2A 3 + A 4 





V f^64 





1 \- 
2 A 54 














2A23 



As one can see, this matrix contains many vanishing elements, and the 8 eigenvalues 
are straightforward to obtain analytically. They read as follows: /_, e 2 , /1, c±, b± 
and pi, where 

Pl = 2(A 3 + A 4 ) - ^A 5 - ^A 6 (3.13) 

As one can see, these additional channels lead only to one extra eigenvalue, p±, 
although we shall see that this eigenvalue plays an important role in constraining 
M H ± and M A . 
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4. Numerical results and discussion 



In this section we present our results for the unitarity constraints on the Higgs 
masses in the THDM. All the eigenvalues are constrained as follows: 

\a±\ , \b±\ , \c±\ , |d±| , |/±| , |ei )2 | , |/i, 2 | , \pi\ < 8vr (4.1) 

In our numerical illustrations, we will consider the special case of A5 = and also 
the general case of A 5 7^ 0. 



Case of A5 = 



For A5 = our potential is identical to those considered in ||14| , [L5| . We improve 
those analyses on two accounts: 

(i) We have considered extra scattering channels which leads to one more eigen- 
value constraint, p±, as explained in Section. 3. 

(ii) When finding the allowed parameter space of Higgs masses we simultaneously 
impose all the eigenvalue constraints. In [|l|| only the condition \a + \ < 8tt was 
applied when deriving mass bounds. 

In order to obtain the upper bounds on the Higgs masses allowed by the 
unitarity constraints we vary all the Higgs masses and mixing angles randomly over 
a very large parameter space. As is reported in [21|, we confirm the result of JK 



which states that a + is comfortably the strongest individual eigenvalue constraint. 
However, the other eigenvalues impose important constraints on and M H ±. If 
only \a + \ < 8tt is imposed we can reproduce the upper bounds on the Higgs masses 
given in [O], in particular their main result of Mh < 410 GeV. When all eigenvalue 
bounds are applied simultaneously we find improved bounds on the Higgs masses, 
particularly for Ma and M H ± . We note that the new eigenvalue constraint pi < 8n 
(eq. |3.13j ) plays a crucial role in determining the upper bound on M H ±. For A 5 = 0, 
we find the following upper bounds: 

M H ± < 691 , M A < 695 , M h < 435 , M H < 638 (GeV) 

Note that the bounds given above are obtained for relatively small tan f3 (say 
tan/3 ~ 0.5). For large tan/3 the bound is stronger, although for the case of A , H° 
and H ± the tan/3 dependence is rather gentle. Of particular interest is the tan/3 
dependence of the bound on Mh which will covered below. 

We have to stress here that our results in the case where the discrete symmetry is 
exact agree with the results obtained using triviality approach 



General case of A5 7^ 

We now consider A5 7^ which corresponds to the inclusion of the term which 
softly breaks the discrete symmetry. Such a term was neglected in the analyses of 
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Figure 1: Maximum Mh in GeV as a function of tan/3 for various values of A5. 



[fL4] , |15| , and from perturbative constraints may take values |As| < 87r p6| . In the 
graphs which follow we do not impose the perturbative requirement |Aj| < 8tt for the 
remaining A». Imposing this condition only leads to minor changes in the numerical 
results which will be commented on when necessary. 

We plot in Fig.l the maximum value of Mh against tan/3 for increasing values 
of A 5 , imposing all the eigenvalue constraints simultaneously as done in Section 3.1. 
For the case of A5 = one finds a strong correlation, with larger tan (3 requiring 
smaller M^. For example, tan/3 > 7 corresponds to Mh < 100 GeV, which is 
the mass range already being probed by LEPII. However, h° in the THDM with 
Mh < 100 GeV is not guaranteed to be found at LEPII due to the suppression 
factor of sin 2 (/3 — a) for the main production process e + e~ — > h°Z. For the case of 
A 5 = we find that values of tan /3 > 20 are strongly disfavoured since they easily 
violate one of the unitarity constraints. If A5 7^ 0, Fig.l shows that for a given tan j3, 
the action of increasing A5 allows larger maximum values of Mh- For A5 = 15 one 
finds a horizontal line at Mh ~ 670 GeV, showing that the upper bound has been 
increased for all values of tan /3. In addition, large (> 30) values of tan /3 are allowed 
if A 5 7^ 0, in contrast to the case of A 5 = 0. For example, A 5 = 1 comfortably permits 



values of tan/3 = 60. However, as pointed out in |27[ perturbative constraints on 
the Aj also restrict the allowed values of tan f3 in the THDM. Using the condition in 
[ 2^| which requires |Aj| < 8tt, we found that tan/3 > 30 is strongly disfavoured. 



We note that the relaxation of the strong correlation between Mh and tan/? 
with A 5 7^ would in principle allow the possibility of distinguishing between the 



11 



discrete symmetry conserving and violating potentials. If h° is discovered and the 
measured values of and tan f5 lie outside the rather constrained region for A5 = 0, 
this would signify A5 7^ and thus a soft breaking of the discrete symmetry. 

5. Conclusions 

We have derived upper limits on the masses of the Higgs bosons both in the Minimal 
Standard Model (MSM) and also in the general Two Higgs Doublet Model (THDM) 
by requiring that unitarity is not violated in a variety of scattering processes. 

In the MSM we have revived the unitarity constraints and shown that the 
inclusion of all the coupled scattering channels in the analysis can make the upper 
limit on the Higgs boson stronger than in the non-coupled case. The same analysis 
has been done for the THDM, where we first considered the THDM scalar potential 
which is invariant under a discrete symmetry transformation and improved previous 
studies by including the complete set of scattering channels. Stronger constraints 
on the Higgs masses were derived. Of particular interest is the tan/3 dependence 
of the upper bound on M h , with larger tan/5 requiring a lighter h° e.g. tan/3 > 7 
implies M h < 100 GeV. We then showed that the presence of the discrete symmetry 
breaking term parametrized by A5 may significantly weaken the upper bounds on the 
masses. In particular, the aforementioned correlation between tan (3 and maximum 
Mh is relaxed. It was suggested that a measurement of tan/3 and may allow 
discrimination between the two potentials. 

Acknowledgements 

The work reported here concerning the THDM has been done in an enjoyable collab- 
oration with Andrew Akeroyd and ElMokhtar Naimi, who are gratefully acknowl- 
edged. 



12 



Appendix A. 

A.l Quartic scalar interactions in the MSM 

The pure quartic scalar interactions in the MSM can be found from the scalar 
potential V = -/i 2 $+$ + A($ + $) where $ = ( , . , r- ] , G ± and G° 

\y V + [(p + iG J/ v 2 J 

are the goldstone boson and 0° is the physical Higgs. Keeping only the quartic terms 
in V we get the following Feynman rules: 

M 2 Q M 2 M 2 

G + G-G + G- = -^- , G + G-G°G° = ?^f , G+CT0O0O *° 



2 tr 2v z 2 v z 

where v 2 is fixed by the electroweak scale by v 2 = (2y / 2Gi?) _1 , M 2 = Av 2 X and 

v 2 = 4 

w 4A- 

A. 2 Quartic scalar interactions in the THDM 



We start from the potential in which the Higgs doublets are expressed in terms 
of the non-physical fields: 



•i i " 



7#< 



IZi 



Expanding this potential and keeping only the quartic terms leads to the fol- 
lowing feynman rules: 



w'lwiwfwi = 4(Ai + A 3 ) , 
Wiwfw2 w£ = 2A 3 + |(A 5 + A 6 ) 
wtw2W^w 2 = 4(A 2 + A 3 ) , 
h 2 h 2 h2h2 = 6(A 2 + A 3 ) , 
hxhiw^wf = 2(Ai + A 3 ) , 
hih 2 wtw2 = (A 8 - A 4 )/2 , 
hih 2 w^w£ = (A 5 - A 4 )/2 , 
z 2 h 1 wf w 2 = |(A 6 - A 4 ) , 



wfwf w 2 w 2 = (A 5 - A 6 ) 
w^wlw^w^ = (A 5 - A 6 ) 

hxhihihi = 6(Ai + A 3 ) 
h 1 hih 2 h 2 = 2A 3 + A 5 

h 2 h 2 Wiwf = 2A 3 + A 4 
h 2 h 2 w 2 w 2 = 2(A 2 + A 3 ) 

hihiW 2 W2 = 2A 3 + A 4 

z 2 z 2 w1[Wi = 2A 3 + A 4 
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z 2 h 1 w±wt = -\ (A 6 - A 4 ) 
z 2 z 2 h 2 h 2 = 2(A 2 + A 3 ) 
z x z x z x z x = 6(Ai + A 3 ) 
z\Z\h 2 h 2 = 2A3 + Xq , 
h\h\z 2 z 2 = 2A3 + Xq , 
z x h 2 w^wi = |(A 6 - A 4 ) 
^i^itu^wj = 2A 3 + A 4 
ziz 2 wfw 2 = (A 5 - A 4 )/2 



-22-22^2"^^ = 2(A 2 + A 3 ) 

-22-22-22-22 = 6(A 2 + A3) 
Z 2 Z 2 Z\Z\ = 2A 3 + A 5 
h 1 h 1 ziz 1 = 2(Ai + A 3 ) 
z x h 2 wtw 2 = -\ (A 6 - A 4 ) 
zxZiwfwi = 2(Ai + A 3 ) 
z 1 z 2 h 1 h 2 = (A 5 - A 6 )/2 
ziz 2 w±w 2 = (A 5 - A 4 )/2 
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